Bose-Einstein Condensation of 84 Sr 
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We report Bose-Einstein condensation of 84 Sr in an optical dipole trap. Efficient laser cooling 
on the narrow intercombination line and an ideal s-wave scattering length allow creation of large 
condensates (No ~ 3 x 10 5 ) even though the natural abundance of this isotope is only 0.6%. Con- 
densation is heralded by the emergence of a low-velocity component in time-of-flight images. 



The study of quantum degenerate gases continues to 
be at the forefront of research in atomic and condensed 
matter physics nearly 15 years after the first observation 
of Bose-Einstein condensation [H, 0, Q . Current areas of 
focus include the behavior of quantum fluids in optical 
lattices Q, effects of dimensionality [H[ and disorder 0, 
0], exploration of the BEC-BCS crossover regime and 
the pursuit of quantum degenerate molecular systems Q . 

Quantum degeneracy in alkaline earth metal atoms and 
atoms with similar electron structure has become another 
area of intense activity. These systems have been the 
subject of many recent theoretical proposals for quan- 
tum computing in optical lattices [id . 
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ation of novel quantum fluids [Hj]. They possess inter- 
esting and useful collisional properties, such as a wealth 
of isotopes that allow mass-tuning of interactions 3, EH 
and creation of various quantum mixtures [l6[ • Low- loss 
optical Feshbach resonances in these atoms 

promise new opportunities because they enable changing 
the atomic scattering lengths on small spatial and tem- 
poral scales. This can lead to creation of matter- wave 
solitons in two dimensions [2(j and random nonlinear in- 
teractions in quantum fluids [2lj]. Many of these ideas 
take advantage of the existence of long-lived metastable 
triplet states and associated narrow optical transitions, 
which are also the basis for recent spectacular advances 
in optical frequency metrology [221 ] . 

Here we report the observation of Bose-Einstein con- 
densation of 84 Sr. To date, quantum degeneracy with al- 
kaline earth metal atoms and similar elements has been 
reported in a collection of ytterbium isotope s |16L l23l . 124 
[25j . and very recently in 40 Ca [26] and 84 Sr [27[, demon- 
strating increasing interest in these systems. Starting 
with a natural- abundance strontium source, our success 
in forming a relatively large Bose-Einstein condensate of 
84 Sr, which has a natural abundance of 0.6%, demon- 
strates the power of laser cooling in strontium using the 
(5s 2 ) 1 So-(5s5p) 3 Pi narrow intercombination line at 689 
nm [28[. It also reflects the near-ideal scattering prop- 
erties of this isotope, which has an s- wave elastic scat- 
tering length of a = 122.7(3) ao [l5j], where ao = 0.53 A. 
The ground state scattering properties of all Sr isotopes 
are well-characterized from one- {2jJ [lo|, HH and two- 
photon [15[ photoassociative spectroscopy, fourier trans- 
form molecular spectroscopy [32[ , and experiments trap- 
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FIG. 1: (color online) Timing diagram for principal exper- 
imental components. The traces are offset for clarity and 
heights schematically indicate variation in power. Note the 
break in the time axis. 



ping various isotopes in optical traps 33, 

Our apparatus and experimental sequence for laser 
cooling atoms and loa ding th em into our optical dipole 
trap were described in [1511351. l3(j ] , and a timing diagram 
is shown in Fig.[H Several techniques are used to enhance 
initial number and subsequent transfers to overcome the 
low abundance of 84 Sr. We form an atomic beam and 
utilize 2-dimensional (2D) collimation and Zeeman cool- 
ing to load a magneto-optical trap (MOT), all using the 
(bs^So-^sbpyP! transition at 461 nm. The MOT con- 
sists of three retro-reflected beams, each with peak inten- 
sity of 3.6 mW/cm 2 and red-detuned from atomic reso- 
nance by 60 MHz. 

One in 10 5 (bsbp^Pi atoms decays to the (5s4a 7 ) 1 D2 
level and then to the (5s5p) 3 Pi and (5s5p) 3 P2 states. 
(5s5p) 3 Pi atoms return to the ground state and are 
recaptured in the MOT, but a fraction of (5s5p) 3 P2 
atoms, with a lifetime of 9min [37(, are trapped in the 
quadrupole magnetic field of the MOT [H, El We 
can accumulate 2.5 x 10 T atoms in an axial magnetic field 
gradient of 60 G/cm with 30 s of loading. (5s5p) 3 P2 
atoms are then returned to the 461 nm cycling transi- 
tion by repumping for 35 ms via the (5s5p) 3 P2-(5,s4a 7 ) 3 D2 
transition at 3012 nm [36]. We subsequently reduce the 
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461 nm beam intensities to about 0.36 mW/cm 2 for 6.5 
ms to reduce the sample temperature to 2mK. 

The 461 nm light is then extinguished, and a sec- 
ond stage of cooling begins using the 689 nm (55 2 ) 1 So- 
(5s5p) 3 Pi intercombination line. This MOT light, con- 
sisting of three retro-reflected beams with 2 cm diam- 
eters, is initially detuned 1 MHz red of resonance and 
broadened by 700 kHz (peak-to-peak dither amplitude) 
to enhance our 689 nm MOT capture rate. The peak 
intensity is 0.75 mW/cm 2 per beam, and the magnetic 
quadrupole field gradient is 0.1 G/cm. Over the next 
150 ms, the field gradient is increased to 0.8 G/cm, the 
laser spectrum is reduced to single frequency and de- 
tuned by only ~ 30 kHz, and the power is reduced to 
0.15mW/cm 2 . This results in 1.6 x 10 T atoms at a tem- 
perature of 1 pK and a peak density of ~ 10 12 cm -3 . 

An optical dipole trap (ODT) consisting of two crossed 
beams is then overlapped with the intercombination-line 
MOT for 115 ms with modest power (2.5 W) per beam. 
The ODT is formed by a single beam derived from a 20 W 
multimode, 1.06 pm fiber laser that is recycled through 
the chamber. No attempt is made to change the polar- 
ization of the beam during recycling. The beams cross 
nearly perpendicularly, but the plane of the lasers is in- 
clined by approximately 10.5° from horizontal. This re- 
sults in a trap with equipotentials that are nearly oblate 
spheroids, with the tight axis close to vertical. Each 
beam has a waist of approximately 100 /im in the trap- 
ping region. This is slightly different than other realiza- 
tions of quantum degeneracy in two-electron atoms, in 
which at least one beam had a waist that was signifi- 
cantly smaller [H, 0, 0, [27] . 

Immediately after extinction of the 689 nm light, the 
ODT power is ramped in 20 ms to 10 W to obtain a 
sample of 3 x 10 6 atoms at 5 pK. The trap depth is 31 
/iK, and the peak density at this point is 4 x 10 13 cm -3 , 
which implies an average collision rate of 1000 s -1 . The 
peak phase space density (PSD) is 10 -2 . For diagnostics, 
we record absorption images of samples released from the 
trap after a time of flight varying from 10 to 40 ms, 
and the optical depth profile can be related to the areal 
density, which, for these long delays, also provides the 
velocity distribution. 

Figure [2] shows the number, temperature, and phase 
space density evolution for a typical forced evaporation 
trajectory. We ramp down the power in the lasers ac- 
cording to P = Po/(l + t/ T )^ 5 with time denoted by t, 
(3 = 1.5, and r = 1.5 s. This trajectory is designed to 
approximately maintain a constant ratio of trap depth 
to temperature (77) during evaporation [39[ . The lifetime 
of atoms in the ODT is 30 s, limited primarily by back- 
ground gas collisions. This allows efficient evaporation 
at a measured 77 = 8 ± 1 for most of the trajectory. 

After 3 s of evaporation to a power of 2W, 1 x 10 6 
atoms remain at a temperature of 0.4 /iK. For this trap, 
our measured trap oscillation frequencies [40[ yield a trap 
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FIG. 2: (color online) Number, temperature and phase space 
density along a typical evaporation trajectory, from 200 ms 
after the completion of loading of the optical dipole trap and 
commencement of evaporation, until quantum degeneracy is 
reached 3 s later. The phase space density increases by a 
factor of 100 for a reduction of atom number by a factor of 3. 



depth, including the effect of gravity, of 3.5 ± 0.3 /iK and 
a mean trap frequency of / = (fxfyfz) 1 ^ 3 — 84 ± 5 Hz. 
These values indicate that the sample is at the critical 
transition temperature for a harmonic trap [411 ]. 



1 KuJN 1 / 3 



(1) 



at which the PSD is 1.2, where N is the number of 
trapped atoms, ks is Boltzmann's constant, Ti is Planck's 
constant, cJ is 27r/, and £ is the Riemann zeta function. 

The evaporation is quite efficient. We lose a factor 
of 3 in the number of atoms from initiation of forced 
evaporation to the onset of degeneracy, while the phase 
space density increases by about a factor of 100. The 
predicted scaling pj pi = (Ni/N) 11 ~ 4 [39], where pi and 
Ni are the initial phase space density and number, and 
7/ = r] + (rj — 5) / (77 — 4), implies that 77 = 7.5, in excellent 
agreement with our estimate from knowledge of the opti- 
cal trap parameters and measured sample temperature. 

Figure [3] shows false color 2D rendering and ID slices 
through the time-of-flight images recorded after 38 ms of 
expansion for various points along the evaporation tra- 
jectory. At 2.7 s of evaporation, the distribution is fit well 
by a Boltzmann distribution, but at 3 s, the presence of a 
Bose-Einstein condensate is indicated by the emergence 
of a narrow peak at low velocity. Further evaporation to a 
beam power of 1.3 W at 4.35 s and trap depth of 600 nK 
yields a condensate with negligible discernible thermal 
fraction. 

The areal density of expanded condensates is fit with 
the functional form 421 
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FIG. 3: (color online) Appearance of Bose-Einstein conden- 
sation in absorption images (left) and areal density profiles 
(right). Data corresponds to 35 ms of free expansion after 
indicated evaporation times (t). (left) Images are 1.8 mm 
per side, and have the same time stamp as the density pro- 
files, (right) The areal density profiles are from a vertical cut 
through the center of the atom cloud, and temperatures are 
extracted from 2D gaussian fits to the thermal component. 
At 3.0 and 3.6 s, a bimodal distribution indicative of Bose- 
Einstein condensation becomes increasingly clear, and then a 
pure condensate is shown at 4.35 s. 



where is the Heaviside function, and R x and R y are 
the condensate radii. At 600 nK trap depth, fitting the 
areal density yields typical values of 7V = 3.0 x 10 5 for the 
largest condensates. The trap mean oscillation frequency 
at this point is / = 70 Hz and harmonic oscillator length 
aho = (h/(MZJ)) 1 ^ 2 = 1.3 /im. The trap harmonic oscil- 
lator energy scale is TuU/kB = 3.3 nK and the chemical 
potential, /i = ^(^) 2/5 = 90 nK. (N a)/a ho = 1500 
is large, implying that the condensate is in the strong 
interaction regime and should be well described by a 
Thomas-Fermi density distribution, confirming that Eq. 
[2] is the appropriate description. 

By releasing a pure condensate from the trap and vary- 
ing the time of flight before recording an absorption im- 
age, we are able to characterize the expansion of the 
many-body wavefunction which is sensitive to the chem- 
ical potential and confinement of the trapped sample. 
Figure [H shows the evolution of condensate radii as a 
function of expansion time obtained by fitting the opti- 
cal depth to Eq. O The fitting coordinate system has 
been rotated to align with the principle axes of the con- 
densate. These samples contain No = 2.3 x 10 5 atoms 
and were obtained after 4.5 s of evaporation to an esti- 
mated trap depth of 350 nK and trap frequencies of 90 Hz 
close to vertical and 60 Hz along the perpendicular axes. 
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FIG. 4: (color online) Condensate radius after release from 
the trap. The expansion of the condensate is fit by a model 
described in the text. Deviations at early times most likely 
reflect limitations of the optical imaging system and large 
optical depths of the condensate. 



The radii show an inversion, as expected for an expand- 
ing condensate, but the effect is small because the trap 
is not highly asymmetric. 

We fit the data using a numerical solution of the 
Castin-Dum model [43[ for the free expansion of a 
Thomas- Fermi wavefunction (42[. The radii evolve ac- 
cording to 



Riit) = Ri(0)bi(t) 



2/i 



m(27r/i) : 



Mt) 



(3) 



for i = x,y, z, where the scaling parameters obey bi(t) = 
(2tt fi) 2 / (bib x byb z ). Using the well-known value of the 
scattering length [l5| and our estimated trap parame- 
ters obtained from beam profiling and trap oscillation 
measurements, we obtain the dashed line in Fig. [4j De- 
creasing the frequencies by about 5% yields an improved 
fit (solid line), which is within our uncertainty in trap 
parameters. 

If the trap is held constant after 4.5 s of evaporation, 
instead of releasing the condensate from the trap, the 
lifetime of the condensate can be measured. Figure [5] 
shows a fit of the decay to an exponential, yielding a 
lifetime of 8 ± 1 s. No discernable thermal fraction was 
evident in these samples. 

The large number of condensate atoms and long life- 
time indicate favorable conditions for many experiments. 
Of particular interest in the near future are the achieve- 
ment of quantum degeneracy with other strontium iso- 
topes and mixtures and the application of an optical 
Feshbach resonance to a quantum degenerate sample to 
modify the scattering length on small spatial and tempo- 
ral scales. 
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FIG. 5: (color online) Number of condensate atoms (filled 
circles) held in the optical dipole trap at constant depth versus 
time after the evaporation (see text). The solid line is an 
exponential fit of the data, yielding a condensate lifetime of 
8±1 s. 
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